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1 Introduction 



We consider an exciton, i.e. an electron-hole pair created in a quantum dot (QD) by 
extremely short optical pulse (Borri 2001, Kikkawa 1988). The rapidly created bare 
exciton is further dressing with phonons, which results in arising a composite quasiparticle 
being the coherent mixture of the electron-hole pair and phonons. Its energy is lower than 
the bare exciton energy. In the case of interaction with only optical phonons this composite 
particle can be called exciton-polaron in analogy to the Frolich electron-polaron. In 
a realistic model much weaker deformation effects responsible for the interaction with 
acoustical phonons arc included. 

The observation of the exciton dressing with phonons is possible when the dipole 
coupling with the electromagnetic field is sufficiently strong to preserve a non adiabatic 
(i.e. rapid) exciton creation (cf. 0.2ps laser pulse in Ref. (Borri 2001)) much faster 
than dressing with phonons (of ps scale as we will show below). We consider a weakly 
polar medium for QDs (GaAs) in which the interaction of electrons with longitudinal 
optical (LO) phonons is dominating in terms of the influence on the quasiparticle energy 
shift. Much weaker interaction with gap-less longitudinal acoustic (LA) phonons plays 
an important role in kinematics of the dressing process, being a channel of energy transfer 
from exciton. In order to determine the limit of decoherence time of the optical excitations 
in QDs, this LA channel for dressing exciton with LO/LA phonons must be included. 

We describe and analyse the dressing and the relaxation of exciton problems in detail 
for typical self-assembled strain induced InAs/GaAs QDs (Jacak 1998) within the Green 
function formulation. We model the QD by the parabolic confinement potential. 



2 Phonon dressing description 

In order to investigate the time evolution of the lowest exciton state, we consider the 
Hamiltonian describing a single electron and a hole interacting with phonons 

x(6k„ + 6Lk„)(e""- - e"'-) - ^ Ll-J^k'/^b^, + 6L^J(a„e""- - ^»e"-'"). (1) 

Here ^^^(feks) denote the phonon creation (annihilation) operators, (s = o for LO, s = a 
for LA), Q denotes the frequency of the zero wave vector LO phonons, Ca - the LA phonon 
frequency, ae,h ~ deformation constants for the electron and hole, M - the mass of ions 
in the elementary cell, v - the elementary cell volume, N - the number of cells in crystal, 
? = (l/^oo — l/eo)"^ - the effective dielectric constant. The electron (hole) part of the 
Hamiltonian 

n^{r,) = + \m*M)\\, + \m\{i^\fzl (2) 

(where i — e,h, = x] + yf) contains the vertical and horizontal confinement potentials. 
There are characteristic confinement sizes /j = ./^-r, kz = W^-r {^o ^D- Finding 
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the eigenfunctions $„(re,r/j) of the exciton Hamiltonian 

He. = Heire) + Hhirh) - r, (3) 

we write the Hamiltonian (|l|) in the second quantisation representation in the form of the 
sum of unperturbed and interaction parts 

H = Y^ hEnttian + ^^i^sblsbks + Fsi.ni, ng, k)4^a„2(6ko + &Lko)' (4) 



where 



Fo(ni,n2,k) = -- ( ^1 j $:,(re,r,)(e^'''-'= -e''''-'')<l>„,(re,r,)d3red3r, 
F,K,n2,k) = - f j $:^(re,r,)M^r^ - a,e''''-'^)$„,(r„ r,)d3r,d3r, 



(5) 



The hnear response of the exciton to the electromagnetic field is described by using 
the retarded Green function of the electric current operators 

~5]d„[a„(t) + 4(t)], (6) 

n 

here d„ denotes an effective dipole moment of the exciton (Elliott 1957, Mahan 2000). The 
spectral density of the linear response function is determined by the causal one-particle 
functions (^^^'(t) = — i(T{an(t)a^/}) (where the average is taken with Hamiltonian (^)) 
which we evaluate solving the Dyson equation. We calculate the spectral function A{ijj) = 
-21mGoo{u) = -2ImG'oo(cu + iO+) and its time course A{t) = (27r)-^ A(cu)e-''^*dcj. In 
order to determine the one-particle causal function, we evaluate the components of the 
mass operator using standard method of the solution of the Green function equations of 
motion 

fi^A _ E^ G^^At) - 1^ E i^s(n,ni,k)i?„^„,(k.,t) = 5(t)5„„,, (7) 



/ ni,ks 



where 



i?„^,.(ks,t) = -i(T{a„,(t)[6ks(t) + fclk.W]«i'(0)}), 



with additional condition {^Pksit)) = (&ks(^) + ^Lks (''')) — (Martin 1959, Engelsberg 
1963). Its Fourier transform takes the form of the Dyson equation 

nin2 

Neglecting the off diagonal Green functions, one finds the diagonal elements of the mass 
operator 



h{uj + uji,s) - - S„i„,(u; + cuks) J 
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■ (10) 



Here n{x) = (e^^'^ — 1) ^, and F^^/ ^^^{uji,uj2) denotes the Fourier transform of a vertex 
function 

, _ 6G-\{t) 
Kn'Msit^ t ) = -F,^n,n\k)6^Ut'y ^^^^ 

satisfying Tnn,ksi^^,^i) ~ 1 — h'^dSnn/duj when consider terms of the lowest order in k. 
Let us denote real end imaginary part of the mass operator, respectively, ReSoo(to') = 
A{uj), ImEoo(cu) = -7(^)- 

We may solve the mass operator equation ( p!o|) with the assumption that the dressed 
exciton is a durable composite particle which is undamped (Suna 1964). It corresponds 
to the strictly solvable one level model of (Krummheuer 2002). We expect in this case 
that the exciton Green function has a real pole at frequency corresponding to the exciton- 
polaron energy. In the frequency region in vicinity of this pole, the imaginary part of the 
mass operator can be neglected compared to the real part, which allows us to write the 
self-consistent equation for (assuming A{uj) —A = const) 

-A^4Ew| + + "^.^^'-"^ I (12) 

A/ ^ [Eq- En- huko Eq- En + huko J 

which may be solved exactly {E^ = Eq — A, J^^,(k) = Fs{n,n' ,'k)Fs{n' ,n, —]<.)). The 
energy shift due to the interaction with LA phonons is negligible. One can calculate the 
form factors 

W = ^(-e - a.)^e-'-'^i/-'^'=^^/^ (13) 



(here l = h^ h = he ^ = U = and = K)^ + ^K^)^ 

m = m* + ml) (Davydov 1976). 

However, we intend to describe the dressed exciton as a damped quasiparticle at non 
zero temperature range. When include many exciton energy levels, the damping is a 
consequence of the inclusion of phonon dispersion which is responsible for the interaction 
between the dressed exciton and the phonon subsystem (Davydov 1972). We solve (|10D in 
two iteration steps in the vicinity of huj = Eq, taking in the first approximation the vertex 
functions F„ ^g{u,uji) = 1, as done in (Moskalenko 1968). Starting from A(o)(co') = 0, 
7(0) (cj) = 0"^ and assuming that A{uj) ~ const, we find in the first iteration step 

A(,)(c.) A(i)(i?oA) = -A' = 



N Eq- En - huJko 

7(1) (cj) = ^ E '^Oo(k) {[1 + n{huJka)]S{huj - Eq - hUl^a) + n{huJl,a) S [tlUJ - Eq + tlUJka)} -(14) 



Here 6„ denotes the n-th level degeneracy rank (6„ = 1,2). The optical phonon contribu- 
tion to 7(^1;), is significant for /icu ~ ± HQ only. We include the LO phonon influence 
on 7(^7), introducing an effective constant to Joo(k), which does not influence the time 
of fast dressing with phonons but it determines the relaxation time at non zero T. This 
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constant will be estimated from the sum rule \ A{t = 0)| = 1 for T = 0. In the second step, 
the relation A(i)(ti;) ^ 7{i)(i^) enables one to neglect the 7{i)(lo') in the denominators of 
(^). It leads to 

HA 

A(2)(^) ^ A(2)(i?o/^) = -A" ^ A^,)iEo/h) + n-'^l^^^^,A' = -A'Z-\ 

7(2)M =7(i)(^ + A'A), (15) 

where En = - A', A,2){Eo/h) = 4 Ek ''"^""^^^ = -A'. Continuing the iteration 
process we would arrive to a self consistent equation for A similar to ([T^), however, after 
infinite number of the approximations. We find 

7(2)(x) = Z-^axh-f^'^'fix) {-n{-x)e{-x) + [1 + n{x)] e{x)} (16) 

where x = hu — Eq + A' , f{x) = J2'^=o ^ ^^(2n+i)n] ' ^ denotes the vertex function 
renormalisation factor {Z = |1 — fi^^dA/dLu\^^). We write the spectral function in the 
form 

where t~^{lj) = Z^iuS) (Langreth 1964, Mahan 1966). Since at zero temperature, the 
renormalisation factor may be calculated directly 

^ ^ N^{E,-En-n^^:)^' ^'^^ 

the sum rule \A{t = 0)| = 1 enables us to estimate the effective damping constant a. At 
higher temperatures, when the phonon dispersion effects are strong we evaluate Z using 
the same condition {\A{t = 0)| = 1). 



3 Results and conclusions 

We have evaluated the interaction form factors i (k) taking the function of the non- 
interacting electron and hole system as the first approximation for the exciton wave func- 
tion ^nive^Th) = ^i';)„^(re)^f,?(re), wherc 

^« (r) = V^?LK,^)0«(^) = y-^j e f e'---^e ^, 



''IZ 



iV« = i ] (19) 

7rV2/. |m|)! )' ^ ' 

Since cug <^ tu*, we consider only the first level for the quantisation in z-axis direction. A 
fact that an effective mass of a hole is large in comparison to the electron effective mass, 
enables one to neglect the hole levels different than the lowest one (of rih = 0, nih = 0). 
Defining the functions 

<'i"(k)= / / ^it(r^,y.)e^^-'^--^(^)V^«(r^,¥.)dr^dy. 



JO 



" (2)e^'=^^0« {z)dz = lf^,^„ (A;x)e-2('=-'-/2)' , (20) 

oo 



where the index n replaces both ne,me quantum numbers, we write the form factors in 
the form 




2 L nn ^ ' nn ^ ' 





hk , 



(21) 



We perform the calculations of the energy shift including ten lowest exciton energy 
levels, for the two confinement sizes 1^ — Gnmfs Ih and le — Snmf^i Ih. The following 
material parameters suitable to the InAs/GaAs QD have been used m* = 0.067mo, mj^ = 
0.38mo, eo = 12.9, = 10.9, hVl = 36.4meV, = 4.8 ■ lO^m/s. The integral form 
factors of the exciton-LO phonon interaction Jq^ = {h / N) J2k JonO^) multiplied by the 
level degeneracy rank corresponding to the exciton energies En are presented in the table. 

The inclusion of many exciton levels influences the dressed exciton energy shift. It is 
responsible for its substantial increase (from 0.5mcV found in one level approximation to 
3.6meV when include ten levels for /g = 3nm, and from 0.35meV to 2.9meV for /g = 6nm). 
The other material constants Ue = 6.7eV, Ue = 2.7eV, p = 5.36g/cm^ (Adachi 1985) are 
useful for estimations. 

The spectral density and its inverse Fourier transform calculated including many exci- 
ton levels and including realistic asymmetry of the QD {Ize/h = 1/3) are plotted in Figs 
1,2. Generally, the LA channel of dressing gives for the typical QDs the picosecond scale 
of dressing. Inclusion of the LO channel does not modify significantly the overall LO 
and LA dressing kinetics in comparison to LA channel solely (Takagahara 1999). Thus 
the LA channel of dressing gives the limit for adiabatic creation of exciton-polaron in 
InAs/GaAs QD. The slope of the time courses of spectral functions is related to damping 
of the dressed exciton at non zero temperature. The relaxation time is strongly influenced 
by the size of QD. 

Let us emphasise that our results for time dependence of the response functions (Fig. 
1,2) coincide with the experimental results by Borri (2001). 



The work was supported by E.C.Project IST-1999-11311 (SQID) and by KBN Project 
PB 2 P03B 055 18. 
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Figure captions: 

Fig. 1. Time course of the response function and spectral intensity for le = 3nm, 

Iz/le = 1/3 

Fig. 2. Time course of the response function and spectral intensity for Ig — 6nm, 

h/le = 1/3 

Table caption: 

Exciton energy levels and corresponding integral form factors Jq„ multiplied by n-th 
level degeneracy rank 6„. 
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